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Abstract. Bender-Canfield showed that a plethora of graph counting problems in orientable/non-orientable 
surfaces involve two constants tg and pg for the orientable and the non-orientable case, respectively. T.T.Q. 
Le and the authors recently discovered a hidden relation between the sequence tg and a formal power series 
solution u{z) of the Painleve I equation which, among other things, allows to give exact asymptotic expansion 
of tg to all orders in 1/g for large g. The paper introduces a formal power series solution v{z) of a Riccati 
equation, gives a nonlinear recursion for its coefficients and an exact asymptotic expansion to all orders in g 
for large g, using the theory of Borel transforms. In addition, we conjecture a precise relation between the 
sequence pg and v(z). Our conjecture is motivated by the enumerative aspects of a quartic matrix model 
for real symmetric matrices, and the analytic properties of its double scaling limit. In particular, the matrix 
model provides a computation of the number of rooted quadrangulations in the 2-dimensional projective 
plane. Our conjecture implies analyticity of the 0{N) and Sp{A'^)-types of free energy of an arbitrary closed 
S-manifold in a neighborhood of zero. Finally, we give a matrix model calculation of the Stokes constants, 
pose several problems that can be answered by the Riemann-Hilbert approach, and provide ample numerical 
evidence for our results. 
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1. Introduction 

1.1. Counting rooted maps in orientable surfaces. This paper introduces and studies the asymptotics 
of a sequence of rational numbers and a conjecture relating them to a sequence of constants (pg) that 
appear in a rooted graph counting problem of [BC| . The problem of counting the number of graphs with a 
fixed number of edges that can be embedded in a surface of genus g has a long history. For planar graphs, 
the problem was solved by Tuttc; [Tu]. In their seminal paper |BC| . Bender-Canfield consider the number 
Tg(n) of rooted maps (G, S), that is embeddings of a graph G with n edges in a closed connected orientable 
surface S of genus g, such that every component of 5* \ G is a disk, and such that an edge, an orientation 
and a side of it is chosen. Bender-Canfield gave an inductive computation of the natural number Tg[n) and 
also proved that for fixed g and large n, Tg{n) is asymptotic to 

(1) Tg{n) ^ t,n^(«-i)A" 
where 

(2) 7 = ^ A = 12 

and the constants tg are computable by some complicated non-linear recursion that depends on an auxiliary 
partition; see [BC[ Eqn.4.2] for a non-linear recursion for 0g*^(/, a) that computes tg via [BC[ Eqn.4.1]. The 
recursion of Bender-Canfield, although cumbersome, gives exact answers for tg that do not depend on any 
unknown constants. In particular, Bender-Canfield obtained the first three values 

_ 2 _ 1 _ 7 

- 7^' ~ 24' - 4320V^- 

In the early nineties, it was realized in combinatorial enumeration that several classes of counting problems 
(such as cubic maps, quartic maps) also lead to an asymptotic expansion of the form ([1]); see for example 
[Gn[G2] . 



1.2. Counting ribbon graphs in orientable surfaces. In the physics literature, graphs appear often as 
Feynman diagrams of a perturbative quantum field theory. In the case of gauge theories with gauge group 
U(A'^), the Feynman diagrams are ribbon graphs, i.e., graphs with a cyclic order of flags around each vertex. 
Perturbative gauge theory counts ribbon graphs with a fixed genus, number of edges and prescribed valency, 
and with weight being the inverse of the size of their automorphism group. This is precisely the content 
of matrix models and was discussed extensively in the eighties, see |BIPZ1 IBIZI . One of the results found 
in this period in the matrix model community was that the generating functions counting ribbon graphs 
at fixed genus are analytic functions with a finite radius of convergence, which is the same for all genera. 
In the late eighties matrix model were studied in the so-called double scaling limit |BKl IDS|, IGM| , where 
roughly speaking one considers the generating functions near their singularity and extracts the coefficients 
Ug of the leading poles for different genera. It was realized that these coefficients are universal, i.e., they 
do not depend on the details of the matrix model potential. It was also realized that the ribbon counting 
problem for arbitrary potential gives rise to a generating series in infinitely many variables that satisfies some 
universal (e.g. KdV or KP) hierarchy. In addition, the double scaling limit of the matrix model satisfies a 
non-linear differential equation which in the simplest case is the famous Painleve / equation; see for example 
|DGJ[ IWij for a survey of these developments in the physics literature. 

Two years ago, Goulden- Jackson proved that a similar generating series that counts maps is a solution to 
the KP hierarchy, and in particular deduced a quadratic nonlinear recursion for Tg(n) in terms of Tg'(ri') for 
suitable {g',n') < {g,n); see fCJ^ Eqn.45]. Using the quadratic recursion of Tg{n), in [BGRj Bender-Gao- 
Richmond give a quadratic recursion relation for the constant tg that involves only tgi for g' < g. This allowed 
them to prove that tg/g^? grows exponentially with a numerical non-zero constant [BGRj . Recently, it was 
realized in |GLM] that the quadratic recursion relation for tg is, in disguise, the recursion characterizing a 
formal power series solution to Painleve I. Although this (as well as Goulden- Jackson's work [Gj:) came as 
a surprise to the enumerative combinatorics community, it is hardly a surprise from the physics point of 
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view. The fact that the tg enter into different map counting problems [Gl| is, from the matrix model point 
of view, a manifestation of universality. 

The analytic structure of solutions to Painleve I is well-known, mostly though the Riemann-Hilbert 
approach (see for example [FIKNj ). and allows one to give the full exact asymptotic expansion of tg in 
inverse powers of g; see for example |GLM| App.A]. 

1.3. The case of non-orientable surfaces. The above discussion focused on the counting problems of 
rooted maps and ribbon graphs in orientable surfaces. Although the problems are different, their view from 
a distance (in the double-scaling limit) is the same, described by a universal non-linear differential equation, 
Painleve I. On the combinatorial side, Bender-Canfield had also developed a theory of counting rooted maps 
in non-orientable surfaces. Let Pgin) be the number of n-edged rooted maps on a non-orientable surface of 
type (7, where 

(3) 

Notice that, in the non-orientable case, g can be integer or half-integer. Then, one has the asymptotics |BC] 

(4) P<,(n)^p<,n^(s-i)A", g > 0. 

This defines a sequence of constants Pg. On the physics side, the Feynman diagrams for 0{N) or Sp(iV) 
theories are ribbon graphs with crosscaps that are embedded in non-orientable surfaces. This is discussed 
in detail in Section [3l If one believes in a matching between the counting problems of combinatorics and 
matrix models, the constants Pg ought to be able to be computed and asymptotically analyzed by the double 
scaling limit of an O(A^) or Sp{N) matrix model. The latter is a pair of functions that satisfies a coupled 
system of two nonlinear differential equations, as was explained in detail by Brezin-Neuberger |BN ' and 
Harris-Martinec |HM| . 

Before we get into details, let us mention that the constants Pg are notoriously hard to compute. In |BC| 
Eqn.3.6] Bender-Canfield give a non-linear recursion for (j)g\l,a) that computes Pg via |BCl Eqn.3.3]. As 
in the case of tg, the recursion is rather cumbersome and gives exact answers for pg that do not depend on 
any unknown constants. In particular, Bender-Canfield obtained the first three values [BC|, p. 245] 

(5) ^i/2--^n7i)' ^3/2 - 3f(T74)' 

The purpose of the present paper is to introduce a sequence (w„) that is easy to compute via a quadratic 
nonlinear recursion relation, and whose asymptotic is easy to analyze, and conjecturally agrees with pg when 
g = {n + l)/2. Using the recursion of Bender-Canfield, Gao was able to match our conjecture for pg for the 
first six values of g. Beyond that, nothing is known. The motivation for the sequence (w„) comes from a study 
of the double scaling limit of O(A^) matrix models following Brezin-Neuberger [BN] and Harris-Martinec 
|HM| . This is discussed at leisure in Section [6l 

1.4. The sequence (m„). Since the counting problems of non-orientable surfaces mix with those of the 
orientable ones (studied in |GLM| ). in this section we give a brief review of Appendix A of |GLM| . Consider 
a function u{z) that satisfies the Painleve I differential equation 

(6) u2 _ i^" ^ 
^ ' 6 

Consider the unique formal power series solution to ^ asymptotic to z^^'^ for large z > 

oo 

(7) 4z) = zl/2^^„^-5«/2^ 

It follows that the sequence (u„) satisfies the following recursion relation 

25(n - 1)2 - 1 1 

(8) U„ = — Un-l ~ 2 ^ UkUn-k, Uq = 1. 

fc=l 
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It was observed that the sequences (ug) are (tg) are related by 

(9) tg^-, ^ 



^9- 29-2r(5^) 

To find the asymptotics of (m„) one uses a trans-series solution of the differential equation For a detailed 
discussion, see Section [5] below. In our context, trans-series are mild generalizations of formal power series 
and can be automatically computed much like the sequence itself. The computation involves a finite 
number of unknown Stokes constants (i.e., adiabatic invariants) S. In our case, we have that the sequence 
(un) has an asymptotic expansion of the form 



(10) Un ~ A- 













1 2m\ 



where the so-called instanton value A and the Stokes constant S are given by 
11 A^——, S^~i-=, 

5 ^/'K 

and the fj,i are defined by the recursion relation 
(12) Mi 



1673/ 



fe=0 



with the understanding that Un/2 = if n is odd. To better understand the recursion relation (I12p . and to 
write it in a more compact form that relates to the trans-series of ([6]), consider the generating series 



oo 



(13) u^{z) = z'^'z-'^/'e-^'-"'" 

Then is equivalent to the following linearized version of ([6]) 

(14) < - 12uMi = 

where u{z) is given by Hopefully, the reader will not confuse the first term ui of the sequence («„) with 



the function ui{z) in Equation (jl4p . 

1.5. The sequence («„). In this section u{z) will denote a function that satisfies the differential equation 
([6]). Consider a function v{z) that satisfies the differential equation 

(15) 2v' -v^ + 3u = 0. 

Consider the unique formal power series solution to pS]) asymptotic to z^/* for large z > 

oo 

(16) z;(z) = ^ i;„z-5"/^ 

ri=0 

It follows that the sequence («„) satisfies the following recursion relation 

(17) Vn = ^-3w„/2 + ^ Vn-l + ^ VkVn~k^ , Vq = -VS 

where (u„) is given by ([5]), with the understanding that Un/2 = if n is not even. Our theorem concerns the 
asymptotics of {v„) for large n. 

Theorem 1. The sequence (w„) has an asymptotic expansion of the form 

(18) Vn ~ (A/2)-"r(n) 
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where A is given in (jlip . S' ^ Q is some nonzero Stokes constant, and the sequence {vn) is defined by the 
recursion relation 

4 ^ 

(19) l^n = -— Vn+l-kVk, VQ = 1- 

on ^ — ' 

An even more compact form of the recursion relation (|19p can be given by introducing the generating 
series 

oo 

(20) i;i(z)=zi/4^i.„z-5n/4^ 

n=0 



Then the recursion relation (|T9|) is equivalent to the following linearized form of (|T5|) 
(21) v[ = vvi 

where v{z) is given by (jl6p . 

1.6. Two conjectures. We can now formulate two conjectures which are motivated by Section |6] below. 
Our first conjecture links the sequence (w„) with the sequence pg. 

Conjecture 1. For all n = 0, 1, 2, . . . we have 

This conjecture reproduces the first three values (O obtained in |BC| and predicts for the next few values 
_ 5 _ 1033 _ 3149 _ 1599895 _ 484667 

^ > P^- 367^' " 1024V6r (f ) ' " 442368' ~ 2949127if7fy ' ~ 56G431872V^' 

Conjecture 2. The Stokes constant S" is given by: 

(24) S" = iVe. 

Remark 1.1. Theorem[T]and Conjecture [T] reveal a single Stokes constant S' associated with the asymptotics 
of the sequence (««). A second Stokes constant is needed for the asymptotics of the fc-instanton expansion 
of (vn); see Theorem [5] and Remark |8. II 

2. Future directions 

2.1. Approaches to Conjectures [1] and [2], In this section we discuss some potential approaches to 
Conjectures [1] and m Recall that Goulden- Jackson construct a solution to the KP hierarchy for a generating 
series associated to the U(A'') gauge group. 

Problem 1. Construct a version of the KP hierarchy corresponding to the the generating series of the 0{N) 
gauge group. 

Perhaps the theory of zonal polynomials will play a role analogous to the Schur functions in generalizing 
the work of |GJj and [O]. This integrable hierarchy without doubt will imply, as in ^GJl Eqn.45] a quadratic 
recursion relation for the number of rooted maps Pg (n) in an non-orientable surface of Euler characteristic 
2 — 2g for a half-integer g. 

Problem 2. State and prove a quadratic recursion relation for Pg{n) in terms of Pg' [n') for {g\ n') < {g, n). 

Given this recursion relation, one may obtain a combinatorial proof of Conjecture [T] along the line of 
thought of (BGRj . A solution to problem [2] can be obtained by identifying Pg{n) with an expectation value 
of an 0{N) matrix model and deduce the quadratic relation from the so-called pre-string equation of the 
matrix model; see for example |MSW| . 

Let us point out that the pair of functions (u, v) have a Lax pair and a Riemann-Hilbert problem, as 
was explained in |BN| . Moreover, it is well-known that every solution to the Riemann-Hilbert problem with 
rational jump functions is meromorphic, see |Mij . In addition, the Stokes constants are exactly calculable 
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from a Riemann-Hilbert problem. For numerous instances of this calculation, that includes the case of the 
Painleve equations, see [FIKNj and also (Kaj . 



Problem 3. Give a solution to the Riemann-Hilbert problem of the pair and compute the Stokes 

constant 5" confirming Conjecture [2l 

2.2. Relation to algebraic geometry. Let Mg^„ be the Deligne-Mumford moduli space of Riemann 
surfaces of genus g with n punctures, and let ■0^, i = 1, • • • , n be the two-cohomology class defined by 

(25) il^i = ci[Ci), z = 1, • • • ,rt, 

where Li is the bundle over Mg^„ whose fiber at Eg e Mg_„ is the cotangent space T*Eg at the i-th puncture. 
In his seminal paper [Wi| , Witten explains how to package the enumerative intersection theory of the moduli 
space of curves into a generating function that ought to satisfy the KdV equation and some initial conditions. 
The conjecture was subsequently proven by Kontsevich [K] . Using the Witten-Kontsevich theorem [Wi[ IK] 
together with the results of TZ] it is possible to show that the intersection number 

(26) (4^-^),^ f_ V?A---A^3V3, 5>1, 
is related to the coefficients Un defined in dH) as follows 



^^^^ (35-3)! (55-5)(55-3)"^- 

A proof of Equation (|77|) is given in [IZl Sec. 6] using the ansatz [IZl Eqn.5.22]. A proof of this ansatz is 
given in 'GJV', Thm.3.1] and also in [EYY] . 

Through © one finds an algebro-geometric interpretation of tg in terms of intersection numbers on Mg^n- 
This relation motivates the following problem. 

Problem 4. Give an enumerative algebro-geometric definition of the numbers Vn (or, equivalently, pg). 

2.3. Acknowledgement. The authors wish to thank J. Gao for an independent computation of the first 
values of the constant Pg using the recursion of Bcnder-Canfield and Nick Halmagyi for computing indepen- 
dently the first values of the number of quadrangulations on the projective plane. M.M. would also like to 
thank Bertrand Eynard for useful discussions. 

3. Analyticity of the 0(iV) and Sp(iV) free energy of a closed 3-manifold 

An initial motivation for our work is the problem of the free energy of a closed 3-manifold M. In |GLMj 
it was shown that the U(iV) version of the free energy F^j{t, h) of M has the form 

oo 
3=0 

where r — Nh and Fj^i g{T) G are formal power series analytic in a disk Dm independent of g that 

contains and depends on M. F^j is defined to be the logarithm of the LMO invariant, evaluated under the 
U(A^) weight system. One may also define the 0(iV) (resp. Sp(iV)) free energy of a closed 3-manifold M by 

(28) (r, h) = W,, (log(ZM)), F^f{r, h) = W,,^ {logiZ,,)) 

where Zm is the LMO invariant of M and Wg denotes the weight system of a metrized Lie algebra g; see 
|B-Nj . The 0{N) (and also the Sp(A^)) free energy of a closed 3-manifold M can be written in the form 

C30 

(29) FZ{r,h)^h-^Y.^9FZ^g{T) 

9=0 

where t = Nh and Fg g{T) e (Q[[t]] for aU g. In [GLMl Rem.4.2] it was observed that if pg satisfies an 
asymptotic expansion similar to tg^ then the 0{N) and Sp(A^) free energy of M enjoys the same analyticity 
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property as the U(iV)-free energy. Thus, Conjecture [T] imphes that the O(A^) and Sp(iV) free energy of M 
is analytic in the above sense. 

Let us discuss in detail the weight system for the Lie algebra Oat of the orthogonal group O(A^). This 
explains the appearance of graphs in non-orientable surfaces. Let Cij denote the square matrix of size N 
with ij entry equal to 1 and all other entries zero. A basis for Oat is Nij = e.y — for 1 < i < j < N . Since 
eijCu = Sj^k^ii (where (5a, 6 = (resp. 1) for a = b (resp. a ^ 1)), and Tr(eij) = Si,j, (where Tr(M) denotes 
the trace of a square matrix M), it follows that the Killing form on o^v is given by 

(30) iN,,,Nki)^2iS,kSa-Sji6^k) 
and its inverse (the so-called propagator) is given by 

(31) {N,,,Nki) = ^{S.kSa-Sjid^k). 

This translates to the following diagrammatic way for the o n weight system of a vertex-oriented cubic graph 
(see for example |B-N| and compare also with [GLMI Eqn.18]): 

i 

^ 1 

2 



The resulting graphs are cubic ribbon graphs with cross-caps, and give rise exactly to embedded trivalent 
graphs in non-orientable surfaces. 

In Section [S] we will study the constants Pg using not the adjoint, but the symmetric representation of 
On- To explain this, observe that if denotes the fundamental representation of Ojv, then A^(C^) is the 
adjoint representation of Otv and Sym (C^) can be identified with the set of symmetric matrices of size N, 
with a basis given by 

Mij — Cij + Cji, for i < i < j < N. 
In that case, the symmetric matrix model (j63p of Section [6] when P — 1/2 has propagator given by 

(32) {M,j,Mki) - ^{SjkS^i + SjiS,k) 

and leads to the following diagrammatic way 

i 

^ 1 

2 



4. The method of Borel transform 

The aim of this section is to prove the existence of an effective asymptotic expansion for the coefficients of 
a non-linear Euler-type differential equation using the theory of Borel transforms; see Theorem [5] below. This 
will provide a proof of Theorem [TJ In fact, Theorem [1] requires only a small portion of the theory, since the 
function v of Theorem [T] satisfies a Riccati type equation (jisp that is well-studied in an excellent exposition 
of [SSj and [BSSi Sec. 5]. Nonetheless, the method of Borel transform is rather general and perhaps not as 
widely known. For the benefit of the reader, we will give a short introduction to this beautiful theory in the 
case of the nonresonant Euler-type differential equation 





(33) 



y' = -Ay (3y + g{x,y) 
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where g{x,y) is analytic at (oo,0), g{x,y) — 0{y^ ,x~^), and the nonresonance condition is A 7^ 0. We wiU 
call A and /3 the eigenvalue and the exponent of (|33p . A and /? are precisely the coefficients of the linearized 
equation 

(34) y' = ~Xy - -fi. 

X 

The Euler-type differential equation includes as a special case the Riccati-type differential equation 

(35) y' = -Xy - ^Py + go{x)y + gi{x)y'^ 

where go and gi are analytic at x = cxd and go{x) — 0{x^^). 

It is easy to see that a nonresonant Euler-type differential equation (|33p has a unique formal power series 
solution /o(2;) G C[[l/a;]]. In fact, Equation (|33p has a unique trans-series solution of the form f{x) where 

(36) f{x)=Y,C''Mx), /fe(a;)=x-^'=e-^^-^//„,fc — , ^1,0 = 1 

k=0 n=0 

is obtained by substituting the above expression in (j33p and equating the coefficient of every power of C 
to zero. This leads to a hierarchy of differential equations for fk{x) which is nonlinear for fc = 0, linear 
homogeneous for fc = 1, and linear inhomogeneous for k > 2. For example, we have: 

(37) /{ = -A/i+g,(x,/o)/i 

The above hierarchy of differential equations was first discovered by Ecalle and was studied in [CNPl p. 54], 
[C] , [Sail p. 11-12] and Sa2, Sec. 8]. An excellent presentation for the special case of the Riccati equation was 
given in [BSSi Sec. 5] and [SSj . This special case is of interest to us since v{z) satisfies the Riccati equation 
()15p . It is important to realize that the coefficients fi„^k are automatically computed from ([33]) by means of 
a nonlinear recursion relation that involves fin',k' for {n' , k') < (n, k). 

The next theorem, presumably well-known to the experts but absent from the literature, links the asymp- 
totics of the trans-series coefficients (/in,fe) in terms of the neighboring coefficients (/i;,fe±i) and two Stokes 
constants ^±1. In the physics literature, fk{x) often occurs as the perturbation theory of the fc-instanton 
solution, in which case Theorem [2] states that the asymptotics of the coefficients of the /c-instanton solution 
can be computed by the (k ± l)-instanton solutions and two adiabatic invariants S±i. 

Theorem 2. Consider the differential equation ()33|) with A 7^ and the coefficients {p-n,k) of the unique 
trans-series solution p6p . Then, for every fc = 0, 1, 2, . . . and n large we have an asymptotic expansion 

(38) - A-"+^(fc + l)^T{n - p) Lo,.+i + E , ^'''^'^ -j 

+ i-xr-^^ik - i)|-ir(n + p) U,^, + f I 

with the understanding that fin^k = for k < 0. S±i are two Stokes constants, defined below. 

Proof. We will give the proof in three steps. At first, let us suppose that fc = and let us denote /i„^o — OLn- 
In other words, we have 

00 

(39) /o(x) = 5]a„— eC[[x-i]] 

ri=l 

is the unique formal power series solution of (j33p . The first step provides the existence of an asymptotic 
expansion of (a„). 

Step 1: With the above assumptions, the sequence (a„) has an asymptotic expansion of the form 



(40) A-"+''|i.r(n-/?) I 
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where Si is an unspecified Stokes constant. 

A proof of Step 1 is given in |CKli Lem.1.1]. The proof of Lemma 1.1 of [CKl] starts by writing down 
a nonlinear recursion relation for the rescaled coefhcients &„ — A"a„/n! of a formal power series solution in 
the form 

bn ~ bn-l + Rn{b{), ■ ■ ■ , ^n-l) 

where i?„(&Oi • ■ • i bn-i) is a nonlinear term, a polynomial in the variables bo, ... , 6„_i. Then, one obtains an 
estimate of the nonlinear term of the form i?n(6o, ■ ■ ■ ,bn) ~ 0{\/n?). This implies that 



It follows that {bn) is increasing and bounded above. Thus, the limit of &„ exists. Using this as input, one 
bootstraps and obtains, order by order in powers of 1/n an asymptotic expansion of (6„). 

The next step identifies the coefficients ni of the asymptotic expansion of (a„) with the coefficients of the 
first trans-series fi{x) of ((33|l . 

Step 2: With the assumption of Step 1, we have 

(41) ^J.n = Mn.l 

for all n. This proves Theorem [2] when fc = 0. 

A proof of Step 2 appears implicitly in pC2l p.l939] and also in p?Pl iBSSl ISSl [Sail [Sa2] . Let us take 
this opportunity to sketch the argument for completeness and also to fix several typographical errors of 
[CK2i p. 1939]. Recall the notion of Borel transform which sends power series in 1/x to power series in p: 

1 

X" ^ r(a) ' 

Consider the differential equation ([55)1 with unique formal power series solution fo{x), and let (j)o{p) denote 
the Borel transform of fo{x) given by 

(42) /o(a:) = ^a„-, ^/.^(p) ^ ^ ^p". 

^ — ' x" ^ — ' n! 

n=0 n=0 

Step 1 implies that <j>o{p) is analytic at p = 0. Consider the punctured plane W — C \ £ where C = AN+ is 
a discrete subset of punctures that lies in a ray A[l, oo). Consider the universal cover W which is identified 
with the set of homotopy classes (rel. boundary) of paths that begin at and lie in M^. Of importance is 
the portion U of W where 

(a) W'^ is the first Riemann sheet of W given by W''^ = C \ A[l, oo). In other words, W'^ is the plane 
minus one cut. 



w 



A 2A 3A4A 

(b) is the second Riemann sheet of W given by all paths in W that cross the cut A[l, oo) at most 
once. 

In |CKli Thm.A.l] it is shown that (po admits analytic continuation as an analytic function in W'^ U and 
(/)o (p) has local expansion around the singularity p = X oi the form 

(43) Ax^a = Si(bi 

where S'l €E C and 

oo 

(«) MP>-p'-'Ef{;^f' 

is the Borel transform of the trans-series solution fi{x) and 

(45) Af,g{p) = lim g{^i + pe'') ~ gi^i + pe~'') 
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denotes the variation of a multivalued analytic germ g near a singularity p = fi; the so-called alien derivative 
in Ecalle's language [Ec| . 

A well-known application of Cauchy's theorem in Borel plane (see jCK2[ p. 1939] and also |CG1|, Sec. 7]) 
now implies the asymptotic expansion (|40p of (a„). Let us sketch the argument, fixing some typographical 
errors from |CK2|, p. 1939], and giving an exact formula, well-known to the physics community. Cauchy's 
theorem for 4>o{p) together with and the analyticity of 0o(p) sd, p — implies that 



I 27ri p"+i ^ 



n 



where 7 is a small circle around p ~ 0. Now enlarge the contour 7 to a contour C U Ti.\ in the first sheet W'^, 
where C is a circle of radius |A| + e, minus an arc, and Ti.\ is a Hankel contour centered at A, for some fixed 




Thus, 



flji+i 



n+1 ^ 



„n+l ' 



nl 27ri J^. 27ri p" 

The first integral can be estimated by 0(n!(|A| -I- e)~"). For the second integral, make a change of variables 
p = A(l + z) and let H denote a Hankel contour centered around 0. Then, we have 



0o( A(l + z)) 



dz = A-" 



Aa'^o(Az) 

+ 



dz 



Now, Equations (gS]) and (|44|) imply that 



1=0 



r(/5 + o 



0+1-1 



and the series on the right hand side, after multiplication by p is analytic at p = 0. A useful Beta- 

integral calculation gives that 

r z-*-' , _ r(7)r(n + i-7) 

Jo + ' 

and therefore 



.7-1 



-dz ■ 



r(n+ 1) 

r(7)r(n-f 1-7) 



(1 + 0((|A| 



')) 



/o (l + z)"+i ^~ r(n-hl) 
Interchanging summation and integration by applying Watson's lemma, (see [Oj) it follows that 



(=0 

Comparing the above with Equation (j40p concludes Step 2. Strictly speaking, the above analysis works only 
when 5R(/3) > — 1. This is a local integrability assumption of the Beta-integral. The asymptotic expansion 
(|40)) remains valid as stated even when 3?(/?) < —1 as follows by first integrating f{x) a sufficient number of 
times, and then applying the analysis. This is exactly what was done in |CK2j and |CG1|, Sec. 7]. 
Let 

00 

denote the Borel transform of fk{x) (given by ([36]) ) for fc = 0, 1, 2, . . . . 
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Step 3: For k — 0, 1, 2, ... , 4>k{p) has analytic continuation as a multivalued analytic function in C \ 
(— fcA + N+A). Moreover, we have 

(46) l\ix<t>k = (fc + l)Si(j)k+i 

for I = 1, —1, —2, —3, . . . , k = 0, 1, 2, 3, ... , k + I = 0, 1, 2, . . . where Si are constants. Moreover [BSS) 
Prop. 5. 4], in the case of the Riccati equation ([35|) . wc have Si ~ for / < —2. Equation ((46|) is the so-called 
Bridge Equation of Ecalle, and appears in [CNPi p. 63] and also in Sal, p. 11-12]. For a detailed discussion, 
see also [BSSl, Prop.5.4] and [SSlIM]. 

Said differently, the Bridge Equation implies that 4>q{p) admits analytic continuation as an analytic 
function in W and the local expansion of every branch of (j)n{p) at each singularity is of the form (|43p for 
a suitable Stokes constant Sj. This is the notion of resurgence coined and studied systematically by Ecalle 
in the eighties. Unfortunately, Ecalle's work remains unpublished, and appears in three Orsay preprints 
that the interested reader is welcomed to read, |Ec| . Fortunately for our purposes, the Bridge equation for 
Euler-type differential equations has been established in print in the above mentioned references. 

The Bridge Equation implies that the nearest nonzero singularity of 4>k{p) appears at p — ±A (resp. 
p = X) when k > 2 (resp. fc = 0, 1), and that the variation of 4>k{p + A) (resp. (pkip — A)) is proportional to 
(j)k+i{p) (resp. 4>k~i{p))- Theorem [5] follows for all k by applying Cauchy's formula to 4ik{p), and deforming 
the contour of integration as in the case of fc = 0. □ 

Let us make some remarks. 

Remark 4.1. As stated, Theorem [5] requires analyticity of the coefficient 

oo 

9{x,y) = ^ .g„(x)y"+^ 

n=-l 

of ([55)1 at (oo,0). In fact, the key Equation (|43p and Theorem [2] remains true when k = and the Borel 
transform 7n(p) of gn{x) has analytic continuation as a multivalued analytic function in C minus a discrete 
set of points, and 7n(p) is analytic in the disk \p\ < |A| -|- e for all n G N. This follows from the proof of the 
resurgence of (t)k{p) in Borel plane; see [CNP| [Sai l [Sa2 l |BSS| ISS] . In particular, Equation gS]) and Theorem 
[2] holds for fc = when the 7„(p) is analytic multivalued in C \ 2Z*A for all n € N, where Z* = Z \ {0}. 
Theorem [2] will be applied in this form to give a proof of Theorem [T] 

Remark 4.2. Starting from a sequence (/in,o) with an asymptotic expansion (|38p . one can consider the finitely 
many sequences (/^n,i) that arise, and repeat this process. This may be continued for ever, arriving at a 
notion of resurgence for sequences (/^n,o) that may not come from differential equations; see for example 
[CGn[CG2l[Gi] . 

Remark 4.3. The constants S±i in Equation psp . if nonzero, can be numerically approximated rigorously 
and efficiently. For examples, see jCCKj IJK[ IBGR] . On the other hand, deciding whether S±i — is 
not algorithmically known. In general the Stokes constants are transcendental invariants of the differential 
equation, sometimes known by the name of adiabatic invariants. 

Remark 4.4. Theorem [2] is valid for systems of first order nonlinear differential equations of rank 1 under 
a nonresonance assumption, see [C] for Steps 1 and 2 and a comment on jSa2| Sec. 13]. In that case the 
nonlinear equation is 

(47) y' = -Ay--By + g{x,y) 

X 

and the linearized equation is 

y' = -Ay - -By 

X 

with eigenvalue matrix A = diag(Ai, . . . , A,.) and exponent matrix B = (/3i, . . . , Pr), where diag(c?i, . . . , dr) 
denotes a diagonal matrix with diagonal entries di, . . . ,dr. The Borel transform (po of the unique formal 
power series solution fo to (|47| is a multivalued analytic function on = C \ £ where C = AiN+ + . . . A,.N+. 
Moreover, 0o has analytic continuation in the principal sheet = C \ Uj^]^Aj[l, oo) which is a plane cut 
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by r rays. The asymptotic expansion of the coefficients of a component of a formal power series solution fo 
of the nonlinear differential equation has an asymptotic expansion of the form 

where J = { j G {1, . . .,r} ||Aj| = min{|A,.| | s = 1, . . . ,r}. 



The analogue of Equation ((38)) does not seem to exist in the literature, especially in the resonant case. 
For some partial results regarding the Painleve I equation, see jCKlj . 

5. Proof of Theorem [T] 

5.1. Theorem [2] implies Theorem [Ij In this section we give a proof of Theorem [1] using Theorem O 
Consider the differential equation where u{z) satisfies ([6]). The change of variables 

(49) u{z)^z^'^{l + U{x)), v{z) = z^'^{-Vi + V{x)), x^?}'^ 

converts ([6]) to the following rank-1 differential equation 

\J'{x) mV(xf 4t/(x) 4 
~x 25 25x2 ~ 25^2 

which can be written as a first order rank-1 differential equation in the form 

(51) U' = -AU- -BU + g(l/a;,U) 

X 

where 



(50) U {x) + — — ^ - ^ = 



25 



and 



g(l/x,U) = 



9m(xf Wix) 4 I 

V 25 25x^ 252;^ / 

A gauge transformation U UG converts (1511) into the normalized resonant rank-1 differential equation 
(53) U' = -AU- iBU-t-g(l/a;,U) 

X 

where 



(^4) A.(^ 'a, B=(2 



1 
1 

2 



and A = 8^3/5 as in ([TT|) . Observe that the differential equation ([53]) is resonant. Nontheless, Theorem [1] 
involves the function v{z) which satisfies a Riccati equation (fT5|) . The substitution (j49p converts psp to the 
following rank 1 differential equation 

(55) n^)-,l^n:.)-^-,^ + ^-^ = o 

5 5a; 5 5x 5 

This is a nonresonant Riccati differential equation of the form ([55]) with eigenvalue A/2 = and exponent 
(3 — Q. Even though the coefficients of ([55)) are not analytic, their Borel transform is a multivalued analytic 
function in C \ Z*y4 and Remark |4 . 1 1 applies . The trans-series solution of (fT5|) is of the form 



(56) v{z)^Y.(^'vi{z) 



1=0 



where ^0(2;) is given in (I16p and 



00 



z;,(z) = zV4e-4-^^^'^^;„_,z-W4 
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and we normalize fo,i — 1- Theorem [T] uses only v\(^z) which satisfies the linearized version 

- -woui = 

of p^ . as stated in pT|) . Theorem[T] follows from Theorem[2]when fc = applied to the Riccati equation ([55 



This Riccati equation does not have analytic coefficients in a neighborhood oi x — oo, however Remark 14.11 
applies. A rigorous numerical calculation of the Stokes constant S' is possible, following standard arguments 
that appear for example in |CCK1 [JK] to conclude that S' ^ 0. This concludes the proof of Theorem [T] □ 

Remark 5.1. One can study in detail the full trans-series solution (j56p . The Vk satisfy the differential 
equations: 

1 

(57) ^fc ~ 2 ^ " ^' 

i=0 

In principle, in writing the equation for the formal trans-series solution to v, we have to include in Equation 
a full formal trans-series solution for u{z), but it is easy to see that the only way to find a solution for 
the Vk is to set this trans-series to zero. From ([57)1 one easily deduces the following recursion relation for 

Vn.l- 

1 StX ^+1 -j^ fc— 1 n+l 

(58) -Wn+Lfe = —\-rV„,k +^'Vn+l^l.kVl + -^^Vi^iVn+l-l,k-i\: fc > 2. 

VS{k- 1) 4 ^ ^=l 1=0 

One finds, for example, 

(59) vo.k = {-l)''-\2Vi)^-\ k>l. 

For equations of the Riccati type, it was shown in |BSS| that the full trans-series expansion can be written 
in terms of three formal power series, as follows. Let us denote 



oo 



(60) Vo{x)^^VnX Vk{x) ^^Vn,kX k > 1. 

Then, there are formal power series v± {x) such that 

(61) Vk^{^lf-^vX-\,'L{l~v+vo), k>l. 
In our case we have 

1 5 1 25 1 3149 1 15995 1 , _g. 

"'^^'^^ ^ ^ 192^/3^ ~ 1152^ ^ 36864^3^ ~ 110592^ ^ ^ 



(62) 



1 1 1 1 1459 1 5429 1 114343 1 , «^ 

V-{x) = 1 h O ix^ ] 

A^/ix 24a;2 11520^3 a;^ 34560 a;^ 138240^3 2;^ ^ ' 



5.2. A brief discussion of the Riemann-Hilbert approach. The Riemann-Hilbert method is an alter- 
native way of proving and computing asymptotic expansions of the form (|48p . For an example relevant to 
the results of the paper, see |Ka l App.A] where Kapaev computes the asymptotic expansion of the sequence 
(u„). The Riemann-Hilbert method uses the fact that all solutions of the Riemann-Hilbert problem (such as 
the function u{z) and conjecturally also v{zy) of isomonodromy are meromorphic functions in the z-plane, 
with prescribed behavior at various sectors. Applying the Cauchy formula and a deformation of the contour 
argument as in Claim 2 above, allows one to deduce the asymptotic expansion of the coefficients (a„) of 
a formal solution; see for example |Ka[ App.A]. In addition, the Riemann-Hilbert approach computes ex- 
actly the corresponding Stokes constants Sj in (|48)) . Whether the Riemann-Hilbert method computes the 
coefficients /i„ in Equation (l4Gp via the recursion p2|) is unknown to us. 

The Riemann-Hilbert approach overlaps with (but is neither a subset or a superset of) the theory of Borel 
transform. On the one hand, the differential equations that are amenable to the Riemann-Hilbert method are 
very closely linked to the notion of integrability, whereas the Borel transform method can deal with generic 
nonlinear differential equations. On the other hand, the Riemann-Hilbert method can deal with variational 
problems that do not come from differential equations. 
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It was recently realized that the method of Borel transform is useful in studying problems of physical origin 
(such as those originating in 3-dimensional Quantum Topology) that do not come from differential/difference 
equations. See for example the results of |CG1[ |CG2] and the survey paper [Gaj . 

Without doubt, the Riemann-Hilbert method is complementary and overlapping with the method of Borel 
transform, and the combination of the two methods can give powerful results which have yet to be witnessed. 

An alternative method to asymptotic expansions, that include exact computation of the Stokes constants 
is available from physics. The part of this theory, relevant to our problem, will be reviewed in Section [71 

6. Matrix models, orthogonal ensembles and non-orientable graphs 

In this section we will give our motivation for Conjecture [1] following the study of the symmetric quartic 
matrix model studied by Brezin-Neuberger [BNj and Harris-Martinec [HM| . 

After reduction to eigenvalues, the so-called partition function of a matrix model can be written as 



(63) ^/^-7OT7^/ndA.|A(A)| 



where 

(64) ^{\i) = \{{h~X,) 

is the Vandermonde determinant of the eigenvalues and V{X) is a polynomial called the potential of the 
matrix model. The index /3 takes the values 1, 1/2, 2 for Hermitian, real symmetric and symplectic matrices, 
respectively. The Gaussian matrix model is obtained for a potential of the form 

(65) V{\) = \\\ 

In order to obtain generating functions of maps, one finds an asymptotic expansion of Z around N = oo. 
One also sets 

(66) gs = -^ 

where t is the so-called 't Hooft parameter and is kept finite in the expansion. This asymptotic expansion is 
called the large N expansion. Since equivalently we are doing the expansion around = 0, we see from (|63p 
that the large N expansion is a generalization of the standard asymptotic expansions of integrals depending 
on parameters. For the case /? = 1/2, which is the one considered in [BN|, IHM] . one finds that F — \ogZ, 
the so-called free energy, can be written in the form 

(67) F{g,,t)^^Foigs,t)+F^ig,,t) 
where Fo{gs^i) and F^{gs^t) have the asymptotic expansions 

oo oo 

(68) F,{g,,t) = ^5s^-2F;(i), F^{g,,t) = ^ gl^~^F^{t), 

9=0 r=l/2 

where in the second sum r takes both integer and half-integer values. The generating functions F°(t) count 
maps on an orientable surface of genus g, while F^(t) count maps on an non-orientable surface with 2r 
crosscaps. These generating functions can be computed in closed form by using elegant methods started 
in [BIPZl IBIZj and which culminated in the algebro-geometric formulation of Eynard and collaborators 
[E1 ICE[ IEO] . It turns out that all these quantities can be calculated by computing residues of meromorphic 
forms on an algebraic curve of the form 



(69) y{x) = M{x)y^a{x), a{x) = {x - a){x - b). 

If V{x) is a polynomial of degree d, M{x) is also a polynomial, of degree d — 2. This curve is called the 
spectral curve of the matrix model, and its detailed form depends on the potential V{z) used in (|63p . Notice 
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that this curve has a branch cut on [a, 6], and the branch points a, 6, as well as M{x) are easy to compute 
once the potential is given (see for example [DGJ[ [M2] for reviews and examples). For a quartic potential 



z2 



4 



(70) V{z) = y + Az 
one has 

(71) b = -a^2a, = ^ (^-1 + VTTiSAi j , 
and 

(72) M{x) = l + 8Xa^ +4,Xx^. 

If we set t — 1 for simplicity, it is easy to see from this description that the generating functions F°{X), 
F^{X) are analytic at A = and that the nearest singularity is at 

(73) A. = -l. 
In terms of the variable y defined by 

(74) =' = .'.-*'°('-^) 

one finds that the generating functions behave as 

(75) 9's'-'F°{X) ^ c°y-5(^-l)/^ g^-'F^W - c^Jy-^C-i)/^, 
for g 7^ 1. For g = I we rather have 

(76) F°iX) ^ 4\ogy. 

The double-scaled generating functions Fl^^{y), F^^{y) are then defined by 

(77) 

r>l/2 

It turns out that these asymptotic expansions can be obtained as solutions to ordinary differential equations. 
For the orientable part, it was shown already in |BK[ !DS, GM] that 

(78) fiy) - -(F^Xiy) 

(also called the double-scaled specific heat) satisfies the Painleve I equation 

(79) f - i/" = y. 

In order to compute the coefficients c° one considers the formal power series for / 

oo 

(80) /(y) = E/«2^"^'""'^^'- 

Tl = 

On the other hand, it was shown in [BNi IHM| that the first derivative 

(81) giy) = -2(^^d" )' 

satisfies the differential equation 

(82) - 6gg' + Ag" - 6gf + 6/' = 0. 

In order to compute the coefficients one considers the formal power series for g 

oo 

5(2/) = E5«2/-^'"-'^/'- 
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The differential equation ([5^ can be integrated out once and gives [BNl Eq.3.25] 

(83) + + c exp (^J' dy' g{y') 

Brezin-Neuberger argue that the appropriate boundary conditions fix c = 0, and it follows that (/, g) satisfy 
the pair of equations 

(84) f^\f"^y, f^^y + lg'. 

Notice that the second equation is a Riccati type equation for g with / known. 

To bring these two equations into the form ^ , we must make the change of variables 

(85) y^2h, f = 2iu. 
The change of variables ([85|) forces the following change of variables 

(86) g^2iv, Tds{z) = 2Fds{y). 

With these change of variables, the differential equations (|84l) for the pair of functions {u, v) become ([6]) and 
psp . In the new variables, the specific heat of the matrix model 

(87) - :F'Uz) = -^(^dJ"(^) - (•^ds)"(^) 
is given by 

(88) u{z) + v'{z) = - E 

ri=0 n=0 

up to an overall factor 1/2. Notice that in the second sum, the index n is related to the half-genus g of |BC[ 
p.244] by 

n = 25 - 1. 

Thus, the total specific heat ([55)) can be written as 



(89) u{z) + v'{z) = Y,[ug~ 



55-3 



^-(53-l)/2 



where g runs now through the natural numbers and the half natural numbers. On the other hand, the total 
heat for this matrix model should behave the same way as the generating series 



g 

of [BCj . The relationship ([9]) between tg and Ug forces us to predict that 

1 5.9-3 _ 1 

which is exactly Conjecture [TJ 

We emphasize that this conjecture is a consequence of the matching between counting problems in com- 
binatorics and matrix models. Although this relation can be rigorously established in certain cases (see 
for example |BDGp . the physics techniques to extract the coefficients tg, pg seem to be more efficient than 
current mathematical techniques, as they lead to simple nonlinear differential equations. 
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7. A PHYSICS DERIVATION OF THE StOKES CONSTANTS 

This section is of independent interest and is included for completeness. It offers an alternative exact 
computation (including the Stokes constant S" of Equation (j24p ) of the asymptotic expansion of sequences 
of enumerative interest, and although it is not rigorous, it is a motivation of Conjecture [2l 

In many quantum field theories, the standard perturbative expansion in powers of the coupling constant 
g has nonperturbative instanton corrections which behave as e~^^/^ ^ where A is the instanton action and k 
is the instanton number. Perturbation theory in powers of 5 at a fixed instanton number k is in principle 
possible, and when successful, it leads to a computation of the coefficients ak,n of series of the form 

C30 

(90) e-'=^/^'g^^afc,„5"- 

n=0 

In the physics community it is believed that, at least in some quantum field theories, the asymptotic expansion 
of {ak,n) for fixed k and large n is exactly computable in terms of the coefficients {ak\n) for fc' nearby to 
k much like Equation ([38]) of Theorem [2l This belief can be mathematically formulated as a resurgence 
conjecture for the power series that appear in perturbative quantum field theory. For an example of this 
principle for a 3-dimensional quantum field theory, see [Gaj . One important aspect of the above belief is 
the ability to compute exactly all constants (including the Stokes constants) in the asymptotic expansions, 
including the Stokes constants, in some class of quantum field theories. This is a reflection of integrability 
of these theories. 

The expansion of gauge theories with gauge group U(A^) is in this respect very similar to the standard 
perturbative expansion in powers of the coupling constant. The expansion also has nonperturbative 
fc-instanton corrections which behave as e"'""'^^ }KNN|, ISh] and have the structure 

00 

(91) e-'^^^^-j:^. 

n=0 

In some cases these series should also have resurgence properties and the corresponding Stokes constants 
should be exactly calculable. 

The field theory relevant to the current paper is a matrix model, where resurgence properties are certainly 
expected. We can study the behavior of the instanton corrections (j90|) in this matrix model near the 
singularity (I73p and extract as before the most singular part of each coefficient afc_„ (which in this case will 
be functions of t,X). One finds that, in this double-scaling limit, the fc-instanton correction ([5U|) . leads to 
the fc-th term in the trans-series expansion of the partition function. Moreover, one can show (at the physics 
level of rigor) that the one-instanton amplitude gives the discontinuity of the first term of the trans-series 
along the Stokes line (see for example [CSj for a detailed exposition), and this makes possible to calculate 
the Stokes constant. In the case of matrix models, this was shown in an important paper by F. David IDaj, 
who calculated explicitly the Stokes constant S in (|11|) by using matrix model techniques. His calculation 
was further clarified in |H-|-| and extended in [MSWj . The connection between instanton calculus in matrix 
models and trans-series of differential equations of the Painleve type is further discussed in [Ml| . In this 
section, we will compute the Stokes constants S and S' by adapting the matrix model technology of [MSW] 
to the case of matrix ensembles with arbitrary values of /3. Needless to say, the computations in this section 
are not rigorous. 

In computing the asymptotic expansion of the partition function (j63p we have expanded around the saddle 
point Ai = for all the eigenvalues i = 1, • • • , N. The first instanton correction corresponds to a saddle-point 
expansion in which one of the eigenvalue integrations is around a nontrivial saddle point X = xq, and can be 
written as 

TV — 1 

(92) /^^^^ ndA.|A(:^,A„...,A.-,)P^e-i^^-^^(^-), 

where the first integral in x is over the saddle-point contour around the nontrivial saddle-point, which we 
have denoted hy x £ X, while the rest of the A'^ — 1 eigenvalues are integrated over the saddle-point contour 
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Xq around the trivial saddle-point. The overall factor of N is due to the fact that there are N choices for x 
among the iV eigenvalues. We have 

4'^(^)-72^^?'(^-l) / d.(det|.l-A.._,PO!;-ie-^^(^) 

(93) 



The notation in this equation is as follows. zj^\N) is the partition function evaluated around the 
trivial saddle-point = 0. A^v-i is the diagonal {N — 1) x (N — 1) matrix given by diag(Ai, • • • , A^r-i). 
(O)^^ is the normalized average of any symmetric polynomial 0{\i) in the eigenvalues A^, computed by a 
saddle-point calculation around the standard saddle-point, 



(94) (O) 



(0) /.,x.nLdA.|A(A)P^0(A)e-^^"^^(^-) 
" /..x„ntidA.|A(A)P^e-^^"^-(^^) ■ 



ueio 

Finally, we have also defined 

(95) fix) = (det \xl - An-iO^Zi e"'^"'^^^. 
The total partition function will be written as a trans-series expansion 

(96) - + 4'' + • ■ ■ 

since, as we will see in a moment, Z^^"* is exponentially small, and it is proportional to the small parameter 
Q--A/gs (where A will be calculated shortly) . It follows that the trans-series expansion of the free energy will 
be given by 

(97) ^ ^(0) _^ ^(1) _^ . . . 
where 

z'>^\n) N 4°^(iV-l) f 

The calculation of this quantity is very similar to the one performed in [MS Wj . and we will just present the 
main intermediate steps. The result involves the connected averages of the matrix model, 

(99) Wh{pi,...,Ph)^p''~'(tr — tr , 

\ pi- An Ph - An / 

where the subscript (c) stands for connected. These averages have an asymptotic gs expansion of the form 

oo 

(100) Whipi, ...,Ph) ^Y.9^s"^'''^WrMpi, ■■■,Ph; t), 

r=0 

where r runs over non-negative integers and half-integers. We will also need the integrated version of these 
averages, 

/Xi j-Xu 
dpi-- - / AphWrMipu--- ,Ph) 
xi—---—Xh—x 

where the integration constant can be fixed by imposing appropriate boundary conditions at x ^ oo (see 
[MSW] for details). Another quantity that enters the computation is the effective potential 

(102) V,six;t)^Vix)~2Ao.,iix;t) 
which satisfies 

/-inoN dVes{x;t) 

(103) = y{x). 

ax 
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We can now compute the integral over f{x) at leading order in gs in terms of these quantities. It is easy to 
see that, due to (|103p . the non-trivial saddle-points are given by the zeros of the moment fimction M(x). 
Let xq be such a zero. We obtain 



(104) 



dxf{x) 



2'Kgs 



exp 



9s 



K=ff(xo) + $(xo) (l + 0(ff.)) 



/ «(2^o) 
where 

(105) <i>{x) ^ (3{2AQ,2ix;t) + dtV,e{x;t)) + Ai/2.i{^;t) 

The quotient appearing in ([98]) can be computed as 



(106) 



r(i + /3) 



exp 



/3 



f3 



9s 2 

Putting everything together, we obtain the contribution of xq to -F'-'^' at leading order in the g^ expansion, 

N r(i + /3) ' 



(107) 

where [MSWj 
(108) 



(2^iV)/3 ^/3(4,r2t)^ V «(^o) 



exp 



/3 



A exp 



<i>(a;o) + ^d^Fo - dtF,/^ 



y{x)dx 



is the instanton action of the matrix model corresponding to the nontrivial saddle-point at Xg. Given a 
matrix model potential, the leading contribution to F'^^^ is given by the sum of the contributions of the 
saddle-points with lowest instanton action (in absolute value). 

All the quantities appearing in (|107p which are needed when (3 = 1 have been already computed explicitly 
in |MSW| . Using the results of [CE we can also compute all the quantities needed for /3 ^ 1, in terms of 
data of the spectral curve ([69]) . Since AI{z) is a polynomial, it can be written as 



d-2 



(109) 

We find for example 

(110) M^i/2,i(p) = (l-r') 
and 

(111) ati^i/2 = -(l-/3)|ilog 



M{z) = cl[{z-z,). 



i=l 



1 2p-a-b 



2^f^ 4ip-a){p-b) 2y^^V P-^ 



/b — a 
V 4 



1/b 

t 



d-2 

logc-h^log -( 



a + b /— — - 

Z^ — + yj(7[z^) 



In the expression for dtFi/2 we have subtracted dtF^r^, which is the r = 1/2 free energy of the Gaussian 
model. 

Let us specialize these results for the quartic matrix model considered in |BN1 IHMj . From (|72|) we see 

that M{x) has two zeros at ±a:o. 



(112) 



Xq 



1 + 8Aa2 
4A 



20 



STAVROS GAROUFALIDIS AND MARCOS MARINO 



Setting /3 = 1/2 one finds by integration of (jllOp that 



A1/2AP) = ^ logb^ - 4a^) - - log[p + Vr^a^ 



(113) 



log 



-log2- - log 



Aa^ - 2x5 - 2a;^ J (1 - 4a2/x§) 



and from (jllip that 
(114) 



1 

dtFi/2 = --log — - log 



i(v/l + 8Aa2 + Vl + 24Aa2) 



Before proceeding with the calculation of (I107p . we present a test of these expressions. Based on general 
arguments relating matrix integrals and enumeration of maps |E2[ IBDGj it follows that 

(115) (trA^)Rp2 = ReSp=o/VFi/2,ib) 

is a generating functional for the number of rooted quadrangulations c„ of the projective plane with n 
vertices, 



(116) 



(trA' 



= t2^c„(-4At)"-i. 



This generating functional can be explicitly obtained from (jll3p in terms of a and xq, 
(117) (tr A^)rp2 =xt-a^- xl{xl + 2^^) Jl - 

\ Xr, 



One finds for the c„ sequence the values 

(118) 5, 38, 331, 3098, 30330, 306276, 3163737, ••• 

The first ones can be tested with the results of |A+| . It should be also possible to obtain the generating 
functional for these numbers from the general combinatorial results of [G3j . 

The zeros (|112p of M{x) give two saddle-points with the same value of A and the same F^^^, therefore 
they both contribute to the instanton amplitude. In order to make contact with the trans-series solution of 
the differential equation, we must analyze the behavior of (|107|) near the singular point ([73]) (we set again 
t = 1). Using the variable z introduced in ([85|) and taking into account the factor of 2 relating F and in 
(l86l). one obtains 



(119) 



3<i7r2 



,-5/6 



exp 



8%/3 



5/4 



: exp 



P-i' ' V2 V 5 
From here, and taking into account the relations between J^{z) and u{z), v{z), we find 

-(1) \iw \ ■ 3~ -5/8 / 8\/3 



(120) 
and 
(121) 



u,{z) = -{T}^lJ'{z) 



2V^ 



v,{z) ^ ~2[T^^lj{z) = iV6zi exp(-^z5/4^ ^ 0{z-^'^)) 



We now recall that this computation computes the discontinuity across the Stokes line, therefore the overall 
constants appearing in these two trans-series solutions are the Stokes constants S*, S' that we introduced 
before. Let us end this section with a problem. 
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Problem 5. Compute rigorously the fc-instanton expansion of a matrix model (with polynomial potential) 
using the Riemann-Hilbert method or the method of Borel transforms. 

8. Numerics 

We now give numerical evidence for Conjecture [2l To do this, we study the sequence 
(122) = 

r(n) 

which according to [1] has the following asymptotic behavior 



oo 



(123) s„~V^-t:, ao = -iS, ai = --S Ai^i, ■ ■ ■ 

fc=0 

In order to test Theorem [1] and Conjecture [2l it is useful to have a precise numerical method to determine 
the numbers a^. The method of Richardson transforms (see for example |BS[ p. 375] for an exposition) is 
very well-suited for sequences of the type (|123p and it is easy to implement. Given a sequence of the form 
p23p . its A^-th Richardson transform is defined by 



k-Q kliN-ky. 

The effect of this transformation is to remove subleading tails in (|123p . The values s^^'' give numerical 
approximations to oq, and these approximations become better as N, n increase. As an example of this 
procedure, we plot the values of the sequences Sn^\ for = 0, 1, 5 and n = 1, • • • , 250. The bottom curve 
corresponds to = 0, the top curve corresponds to = 1, and the intermediate curve is A^ = 5. 




We can see that the Richardson transforms provide a very fast convergence to the expected value 
(125) = 2.44948974278317809819728407471.... 

Numerically, we find that n = 250 and A^ = 20, 30 provide an approximation to the expected value with 28 
(respectively, 30) significant digits, 



(126) 



sg°J = 2.44948974278317809819728407459..., 
sg°) = 2.44948974278317809819728407471... 



We can verify in a similar way the value of oi by considering the auxiliary sequence 
(127) r„=n(i|i-l) 
with the asymptotic behavior 



(128) 



r„ ^ ^i^iA + ©(l/n), n^oo. 
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We plot the values of the sequences rl^\ for iV = 0, 1, 5 and n — 1, - ■ ■ , 250. 




Like before, we have very fast convergence to the expected value 

-viA = — . 
2 5 

„(20) 



(129) I'^iA^-l- 



One finds numerically 
(130) 



^250 = -0.200000000000000000000000001520..., 
r^go^ = -0.200000000000000000000000000002.... 



Remark 8.1. We can also study the asymptotics of the coefficients Vn,k appearing in the trans-series solution 
([5^)1 . This asymptotics is governed by Equation ([55]) . where A = A/2, /3 = 0, and Si — S' . Based on our 
numerical results, we conjecture that 

(131) = -i^. 

Finally, we end with the following expectation of Conjecture [TJ 

1238878081129358302459331398309144842472024202171957968278854904568087551305256373 



Pii = 



10986030082548950321157435333449889551411576832V6r (if^) 
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